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摘 要 本文 基于 非 平衡 分 子 动力 学 方法 ， 模 拟 研究 了 梯形 、 三 角形 和 甜 形 热 脉 冲 在 石墨 烯 中 的 传播 过 程 和 规律 。 三 种 类 型 的 
温度 脉冲 采用 直接 速度 修正 法 ， 通 过 调整 对 控 温 区 原子 施加 的 温度 增 量 来 实现 。 结 果 表明 ， 热 扰动 在 石墨 燃 中 以 有 限 的 速度 传 
播 ， 并 表现 出 明显 的 波动 传递 特征 。 ni 道 波 具有 宏观 动量 ， 波 速 等 于 
石墨 烯 声 子 群 速度 ， 为 声波 ， 另 一 道 波 是 无 宏观 动量 的 热量 的 传播 ， 波 速 是 为 热流。 研究 还 发 现 ， 热 脉冲 形状 
对 热 波 现象 的 产生 以 及 声波 、 热 波 的 传播 速度 没有 影响 ， ER 
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Molecular Dynamics Simulations of Heat Pulse Propagation in Graphene 
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Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China) 


Abstract: Based on nonequilibrium molecular dynamics, heat pulse propagation in graphene is studied by imposing 
trapezoid, triangle or rectangle shaped heat pulses. The three kinds of heat pulses are established by velocity scaling 
methods. It is found that the heat pulse in graphene propagates with finite velocity and excites two waves. One of the 
two waves 1s a Sound wave, which has macroscopic momentum and propagates at the Speed of sound. The other ls a 
thermal wave, whose propagation speed is 1/V3 of the sound velocity. The shape of the heat pulse has no influence 
on the thermal wave occurrence and velocities of the two waves, but it remarkably Influences the peak temperature. 
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0 引言 豫 范 围 ， 导 热 过 程 偏离 傅立叶 导热 定律 ， 出 现 瞬 态 
非 傅立叶 导热 或 热 波 现象 EC。Osman 等 局 模拟 了 碳 

石墨 烯 是 由 碳 原子 紧密 排 布 成 的 六 边 形 二 维 材 纳米 管 在 低温 下 的 瞬 态 导热 过 程 ， 发 现 热 脉冲 先后 

料 ， 自 2004 年 首次 成 功 制备 以 来 ， 石墨 烯 凭借 其 优 。 激发 出 三 个 波 包 ， 它 们 分 别 对 应 纵向 模式 、 卷 曲 模 

异 的 结构 和 物理 性 质 迅 速成 为 科学 研究 和 工程 应 用 式 声学 声 子 和 第 二 声 。 Shiomi 等 四 在 碳 纳米 管 中 模 拟 

的 关注 热点 趾 。 其 中 一 个 重要 的 研究 方向 就 是 利用 其 。” ”观察 到 了 热 波 ， 并 结合 模 态 分 析 判 定 声学 声 子 以 弹 

超 高 的 导热 和 导电 特性 制备 以 石墨 烯 为 基础 材料 的 。”” 道 输 运 方式 传递 ， 光 学 声 子 对 非 傅立叶 导热 有 主要 

微 / 纳 电子 器 件 , 解决 芯片 在 高 集成 度 下 的 散热 问题 ， 贡献 。 

从 而 确保 其 运行 的 可 靠 性 外。 尽管 矶 纳米 管 与 石墨 烯 同属 低 维 碳 素材 料 且 具 
对 于 一 般 情况 下 的 导热 问题 ， 健 立 叶 导热 定律 ” 有 很 多 相似 之 处 ， 但 两 者 在 导热 性 能 方面 仍 存在 差 

有 足够 的 描述 精度 。 但 对 于 微 /纳米 尺度 下 的 超 快速 异 。 目 前 关于 石墨 烯 瞬 态 导热 特性 的 研究 还 很 少见 

导热 问题 ， 其 时 间 和 空间 尺度 接近 或 小 于 材料 的 弛 到 报道 ， 因 此 本 文 以 石墨 烯 为 研究 对 象 ， 采 用 非 平 
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衡 分 子 动力 学 方法 对 石墨 烯 中 热 脉冲 的 传播 过 程 进 
行 了 模拟 ， 并 在 此 基础 上 对 比 总 结 了 不 同类 型 热 脉 
冲 的 传播 特性 和 规律 。 

1 模拟 方法 

本 文采 用 非 平衡 分 子 动力 学 方法 模拟 了 石墨 烯 
的 瞬 态 导热 过 程 ,模拟 系统 如 图 1 所 示 。 系 统 长 99.77 
nm、 宽 2.16 nm, 长 、 宽 方向 均 采 用 周期 性 边界 条 件 。 
碳 碳 原子 间 的 作用 力 采用 Brenner 势能 函数 描述 上 
即 : 

E -> (ns) -bv (| (1) 
其 中 ，E。 是 总 势能 ，ftry) 是 势能 截断 函数 ，VR、VA 
分 别 是 原子 间 的 排斥 势 和 吸引 势 ，b;j 是 多 体 作 用 系 
数 ， 表 示 周 围 原子 对 i、j 原子 作用 势 的 影响 ， 与 原 
子 间 键 角 有 关 。 系 统 原子 相 空 间 的 演化 采用 
Leapfrog-Verlet 格式 算法 ， 积 分 步 长 取 0.5 信 。 将 模 
拟 系 统 沿 x 方向 等 分 为 200 个 区 域 ,根据 能 量 均 分 原 
理 ， 各 区 域 的 瞬时 温度 可 以 通过 统计 平均 原子 动能 
获得 : 


2 :二 .BY? 
T th 0 py 
WD) 7m 2 


其 中 ，T(,D 是 温度 ，m 是 矶 原子 质量 ，N 是 各 区 域 
原子 数 , ks 是 玻 尔 兹 曼 常 数 , Pi 是 第 i 个 原子 的 动量 。 
需要 特别 说 明 的 是 ， 当 系统 处 于 非 平衡 状态 或 存在 
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图 1 中 (a)、(b)、(o) 所 示 的 梯形 、 三 角形 和 和 矩形 脉冲 
的 函数 分 别 为 : 


~ f1000, 0<t<0.5 
? 1-1900t+1950,.0.5<t<1.0， Ga 

ieoot+50， 0x<t<05 
? 1-1900t+1950,0.5<t<1.0” Cy 
T, =1000，0<t<1.0， Go) 


其 中 Tp 是 脉冲 温度 ， 单 位 为 K，t 是 时 间 ， 单 位 是 
ps。 三 种 脉冲 最 高 温度 均 为 1000 KK， 持续 时 间 均 为 
1.0 ps。 以 开始 施加 热 脉 冲 的 时 刻 为 零 时 刻 ， 加 热 停 
止 后 系统 在 绝热 条 件 下 演化 ， 至 4 ps 时 模拟 结束 ， 
此 时 系统 中 心 区 域 的 热 扰 动 已 传播 至 长 度 方向 的 边 
界 。 另 外 ， 为 消除 噪音 对 模拟 结果 的 影响 ， 本 文 改 
变 初始 参数 进行 了 20 组 模拟 ， 由 系 综 平均 获得 系统 
的 温度 及 动量 分 布 。 由 于 模拟 系统 在 长 度 方向 是 对 
称 的 ， 所 以 本 文中 只 讨论 0-50 nm 的 区 域 。 
2 结果 及 讨论 

图 2 以 温度 云图 的 形式 给 出 了 梯形 热 脉 冲 ( 式 3a) 
在 石墨 烯 中 的 传播 过 程 。 这 里 横 、 纵 坐标 分 别 表示 
空间 位 置 和 时 间 ， 颜 色 深浅 示意 温度 的 高 低 。 云 图 
中 热 脉 冲 施加 在 x=0~1.5 nm 和 扩 0~1.0 ps 的 区 域 ， 
在 图 中 以 矩形 虚线 框 表示 ， 下 文中 将 该 区 域 称 为 控 


宏观 运动 时 ，Ttxb 仅 代表 区 域内 原子 的 平均 动能 ， 
但 为 方便 讨论 ， 本 文 仍 将 其 表示 为 温度 。 


nm 


图 1 模拟 系统 示意 图 


Fig.1 Schematic of the simulation System 


初始 0.5 ns 采用 Nose-Hoover 热 浴 使 系统 达到 平 
衡 温度 50 K， 之 后 采用 直接 温度 修正 法 在 系统 中 间 
区 域 (图 1 阴影 所 示 区 域 ) 施 加 持续 时 间 1.0 ps 的 温度 
脉冲 。 有 具体 方法 是 采用 速度 修正 因子 (TyD 1 调整 原 
子 速度 ， 使 中 间 区 域 的 瞬时 温度 达到 指定 脉冲 温度 。 


温 区 。 


Xinm 


鸡 


2 施加 梯形 热 脉 冲 的 系统 温度 云 医 


Fig.2 Temperature profile of the system imposed a trapezoid 


shaped heat pulse 


从 图 2 可 以 看 出 ， 如 前 所 述 ， 对 平衡 系统 温度 
场 的 扰动 是 通过 对 控 温 区 原子 施加 温度 增 量 构造 
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的 。 所 施加 的 热 脉 冲 随 着 时 间 的 推演 由 控 温 区 癌 系 
统 两 端 传 播 ， 图 中 虚线 表示 扰动 传播 的 波 前 位 置 ， 
其 在 3.3 ps 左右 到 达 系 统 边界 , 这 充分 说 明 在 本 文 的 
模拟 中 ， 热 扰动 的 传播 速度 是 有 限 值 。 另 外 ， 在 传 
播 过 程 中 ， 热 脉冲 共 激 发 出 两 道 以 不 同 速度 向 边界 
传递 的 波 , 其 中 第 一 道 波 ( 实 线 1) 速 度 较 快 , 在 3.5 ps 
左右 到 达 系 统 边界 ， 第 二 道 波 ( 实 线 2) 速 度 较 慢 ， 在 
4ps 时 仅 传 播 至 35 nm 附近 。 上 述 扰动 传播 速度 的 有 
限 性 以 及 波动 传递 特征 均 是 瞬 态 非 傅立叶 导热 的 体 
现 。 

傅立叶 导热 定律 导出 的 温度 演化 方程 会 得 出 无 
限 大 的 热 扰 动 传播 速率 ,为 克服 这 个 问题 , CattaneoDl 
和 Vernotte 等 通过 引入 热流 对 时 间 的 导数 项 提出 了 
CV 模型 : 


0q 
qtr = VT, (4) 
其 中 gq 是 热流 密度 ，z 是 弛 了 移 时 间 ，k 是 热 导 率 。 
CV 模型 使 温度 演化 方程 变 为 双 曲 型 波动 方程 ， 从 而 
可 以 得 到 有 限 的 热 扰 动 传播 速度 ， 即 : 


v=v, /V3, (5) 


这 里 vy 是 热 波 速度 ，vg 是 声速 。 

图 3 给 出 了 波峰 在 不 同时 刻 的 位 置 ， 可 以 看 出 
热 脉冲 激发 的 两 道 波 在 传播 过 程 中 波 速 不 变 ， 拟 合 
直线 可 得 w=13.4 km/s、v=9.33 km/s。 由 文献 [10] 可 
知 模拟 结果 中 第 一 道 波 的 速度 与 石墨 烯 声 子 群 速度 
基本 相等 ， 并 且 第 二 道 波 与 第 一 道 波 的 速度 也 基本 
符合 WV3 倍 的 关系 。 


1.8 21 24 27 30 33 3.6 3.9 


3 波峰 位 置 随时 间 的 变化 


Fig 3. Location of wave crests varying with tme 


图 4(a) 给 出 了 二 1.2 ps、2.0 ps 和 2.8 ps 时 石墨 烯 
的 温度 分 布 以 及 2.8 ps 时 的 动量 分 布 。 图 4(b) 给 出 了 
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相同 时 刻 由 去 掉 局 部 宏观 动量 的 原子 速度 计算 的 真 
实 温度 和 动量 分 布 。 

从 图 4 可 以 看 出 ， 热 脉冲 激发 的 两 道行 波 在 传 
播 过 程 中 能 量 逐 渐 衰 减 。 同 时 受热 扰动 影响 的 区 域 
均 存 在 宏观 动量 ， 且 0-30 nm 区 域 与 30-45 nm 区 域 
的 局 部 动量 方向 相反 ， 系 统 整 体 表 现 出 纵向 疏 密 波 
的 特征 。 在 去 掉 局 部 宏观 动量 后 ， 第 一 道 波 完全 消 
失 ， 第 二 道 波 波峰 温度 有 所 衰减 。 以 上 结果 均 表 明 ， 
施加 在 控 温 区 的 热 脉冲 在 传播 过 程 中 激发 的 第 一 道 
波 是 系统 受热 膨胀 所 形成 的 声波 ， 第 二 道 波 是 由 唱 
格 振动 形成 的 温度 波 ， 即 热 波 。 


5 10 15 
xX/nm 


SS 10 '15 -20. :25 
x/nm 


图 4 不 同时 刻 温度 和 动量 分 布 (a) 有 宏观 动量 ;(b) 无 宏观 动 


[= 
时 


Fig 4. Temperature and momentum profiles at different time (a) 
with macroscopic momentum; (b) without macroscopic 


momentum 


图 5(a)、(b) 分 别 为 施加 三 角形 ( 式 3b) 和 和 矩形 ( 式 
3c) 热 脉冲 时 的 温度 云图 。 比较 图 2 和 图 5 可 以 发 现 ， 
三 种 超 快速 热 脉冲 在 纳米 尺度 石墨 烯 中 的 传播 过 程 
均 表 现 出 明显 的 波动 传递 特征 ， 并 且 脉 冲 形状 的 改 
变 对 其 激发 的 行 波 波 速 和 宽度 没有 影响 ， 但 对 峰值 
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温度 影响 显著 。 在 3 ps 时 ， 梯 形 脉冲 对 应 声波 


ChinaXiv 合 作 期 刊 


所 形成 的 纵向 琉 密 波 ， 为 声波 ， 男 一 道 波 无 宏观 动 


的 峰值 温度 为 465 K, 三 角形 为 257 K, 算 形 为 326K。 ” 量 ， 其 速度 与 声速 基本 符合 WV3 倍 的 关系 ， 是 由 无 
在 三 种 脉冲 中 ， 和 矩形 脉冲 能 量 最 高 ， 梯 形 次 之 ， 三 岗 则 蝇 格 振动 所 形成 的 温度 波 ， 为 热流。 通过 比较 
角形 最 小 。 可 见 脉 冲 激发 的 行 波峰 值 能 量 的 大 小 不 ”三 种 热 脉 冲 的 传播 过 程 发 现 ， 热 波 现象 的 产生 以 及 
能 直接 由 热 脉 冲 能 量 的 多 少 决定 ， 也 与 热 脉 冲 在 高 。 行 波 的 波 速 均 不 受热 脉冲 形状 的 影响 ， 但 其 峰值 温 
温 段 的 持续 时 间 有 关 。 高 温 段 持续 时 间 短 且 能 量 高 。 度 与 脉冲 形状 有 关 。 

的 脉冲 更 容易 激发 出 共有 较 高 峰值 能 量 的 行 波 。 
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5 施加 (9) 三 角形 、(b) 和 矩形 热 脉冲 的 系统 温度 云图 
Fig.5 Temperature profiles of the system imposed (a) a triang 
and (b) a rectangle shaped heat pulse 


3 结论 


本 文采 用 非 平 衡 分 子 动力 学 模拟 方法 ， 通 过 


le 


向 


石墨 烯 施加 持续 时 间 为 1 ps 的 梯形 、 三 角形 和 拢 
热 脉冲 ， 研 究 了 热 脉冲 在 石墨 烯 中 的 传播 规律 。 
扰动 的 传播 过 程 通过 统计 平均 时 空域 上 的 原子 动 
温度 描述 。 结 果 表 明 ， 石 墨 烯 中 热 扰动 以 有 限 速 
传播 ， 并 表现 出 明显 的 波动 传递 特征 。 热 扰动 激 


I 


能 
度 
发 


县 ， 


出 两 道 以 不 同 速度 传递 的 波 ， 一 道 波 具 有 宏观 动量 


其 速度 与 声 子 群 速度 基本 相等 ， 是 由 系统 受热 膨 


胀 
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